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Abstract. Some theories of new physics accounting for neutrino masses can give rise to a low-energy non-unitary leptonic 
mixing matrix. It is shown that the CP-asymmetries in the — ► v T channel are an excellent probe of such new physics. In 
adition, we clarify the relationship betweeen our framework and the so called "non-standard neutrino interactions" scenarios: 
the sensitivities explored here apply as well to such constructions, except for extremely fine-tuned cancellations. 
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INTRODUCTION 

Non-zero neutrino masses are evidence for physics beyond 
the Standard Model (SM). The complete theory account- 
ing for them and encompassing the Standard Model should 
be unitary, as mandated by probability conservation. The 
light fermionic fields in the theory may mix with other de- 
grees of freedom. The complete mixing matrices are genet- 
ically unitary. However, the effective 3x3 submatrices de- 
scribing the mixing of the known fields need not to be uni- 
tary. We will assume that the full theory is indeed uni- 
tary, whereas, low-energy non-unitarity may result from 
BSM physics contributing to neutrino propagation, when 
the physical measurements are described solely in terms of 
SM fields 1 1]. Specifically, the tree-level exchange of heavy 
fermions (scalars) will (not) induce low-energy non-unitary 
contributions |2]. 

In Ref. |1] the so-called MUV (minimal unitarity vio- 
lation) was developed and the absolute values of the ele- 
ments of the matrix N were determined. However, no in- 
formation on the size of the phases of the mixing matrix is 
available, neither on the standard "unitary" phases nor on 
the new non-unitary ones, as present oscillation data cor- 
respond mainly to disappearance experiments. It is the pur- 
pose of this work to explore the future sensitivity to the new 
CP-odd phases of the leptonic mixing matrix associated to 
the non-unitarity. Furthermore, a comparison between our 
approach and the results in the "non-standard neutrino in- 
teraction" scenarios 
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FORMALISM 

Let us parameterize the general non-unitary matrix N, 
which relates flavor and mass fields, as the product of an 
hermitian and a unitary matrix, denned by 



with tj ' = tj. The bounds derived in Ref. 1 1] for the modu- 
lus of the elements of NN^ also apply to the elements of 7], 
since AW"'" = (l + rj) 2 ^l+2r; and it follows that 
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at the 90% confidence level. The bound on tj^ t has been 
updated with the latest experimental bound on t — > py | 6]. 
Eq. © shows that the matrix A' is constrained to be unitary, 
within a 10~ 2 accuracy or better. The unitary matrix U in 
Eq. |Q~|| can thus be identified with the usual unitary mixing 
matrix U = Upmns* within the same accuracy. The flavor 
eigenstates can then be conveniently expressed a£] 
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It follows that the neutrino oscillation amplitude, neglect- 
ing terms quadratic in r\ , is given simply by 



< Vp\v a (L)>=A S a M R (L) (l-naa-Tjpp) 



with 
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xpK^J =^ v p \"a 
being the usual oscillation amplitude of the unitary analy- 
sis. 

New CP-violation signals arising from the new phases 
in 7J, require to consider appearance channels, a ^ p. The 
best sensitivities to such phases will be achieved in a regime 
where the first term in Eq. 10 is suppressed. This happens 
at short enough baselines, where the standard appearance 



The handy superscript SM is an abuse of language, to describe 
the flavor eigenstates of the standard unitary analysis. 



amplitudes become vanishingly small while A^E(L) ~ 1. 
The total amplitude is then well approached by 

< v p | v a (L) >=A^p(L) +277*^ + 0(r\A), (6) 

where ff(r\A) only includes appearance amplitudes and 77 
components with flavor indices other than a/3. Then, at 
short enough baselines, each oscillation probability in a 
given flavor channel, P a p, is mostly sensitive to the cor- 
responding T] a p . The other elements of the 77 matrix can be 
safely disregarded in the analyzes below, without imply- 
ing to assume zero values for them. That is, their effect is 
generically subdominant, a fact that has been numerically 
checked for the main contributions. 

For instance, in a two family scenario and within the 
above-described approximation, the oscillation probability 
would read: 

P a/3 =sin 2 (20)sin 2 ^J 

-4|r?«/j|sin5 a/ jsin(20)sin^j +4|r7 a/ j| 2 , (7) 

where A = Am 2 /2E and r] a p = | r] a p |e~' 5 »P . The first term 
in Eq. is the usual oscillation probability when the 
mixing matrix is unitary. The third term is the zero-distance 
effect stemming from the non-orthogonality of the flavor 
eigenstates. Finally, the second term is the CP-violating 
interference between the other two. Notice that, even in two 
families, there is CP-violation due to the non-unitarity. 

In the section below, we will analyze the new sources 
of CP violation in the V/i — > v T channel, since present 
constraints on r\ el i are too strong to allow a signal in 
the v e — * Vfi one (see Eq. l(2j), and v e — > v T has extra 
supressions by small standard parameters such as sin 613 
or A12 LZy. When numerically computing P^ x , the only 
approximation performed will be to neglect all 77 elements 
but r\fx t . They should be indeed subdominant, as illustrated 
by Eq. l|6j. In any case, we have checked this approximation 
numerically. 



SENSITIVITY TO THE NEW CP-ODD 
PHASES: THE v M -> v T CHANNEL 

As suggested by Eqs. Q, the best sensitivities to CP- 
violation will be achieved at short baselines and high en- 
ergies, where the standard term is suppressed by sin 2 (4p). 
We will therefore study a Neutrino Factory beam 1 8] result- 
ing from the decay of 50 GeV muons, to be detected at a 
130 Km baseline, which matches for example the CERN- 
Frejus distance. For these values, sin(^y^) ~ 1.7 • 10~ 2 and 
sin(^) ~ 6- 10~ 4 , where A^ = (mj-mf) /2E. All terms 
in the oscillation probability Eq. (|7) can then be of similar 
order for the channel v)i — > v T , if the 77^ x value is close 
to their experimental limit in Eq. l(2j. In what follows, we 
will assume 2 ■ 10 20 useful decays per year and five years 



running with each polarity, consider a 5 Kt Opera-like de- 
tector and, finally, sensitivities and backgrounds a factor 
5 larger [7] than those used for the v e — > v T channel in 
Ref. 11. 

In Ref Q| the complete expanded expression for/^ T can 
be found. However, as this probability is not suppressed by 
small standard parameters such as sin 613 or A12, the two 
family approximation in Eq. Q, with 8 = 023 and A = A31, 
is very accurate to understand qualitatively the results. That 
equation indicates that the CP-odd interference term is only 
suppressed linearly in |t)ut|. This can indeed be observed 
in the result of the complete numerical computation, Fig.[T] 
which shows the sensitivities to | T7 jti - r | and S !1T obtained. 
The left panel represents two fits to two different input val- 
ues of 1 17/1 t I and 5/i T (depicted by stars). The dashed lines 
correspond to fits done assuming the wrong hierarchy, that 
is the opposite sign for A3 1 to that with which the number of 
events were generated. As expected from Eq. (0, a change 
of sign for the mass difference can be traded by a change of 
sign for S^ t . Nevertheless, this does not spoil the potential 
for the discovery of CP violation, since a non-trivial value 
for |5ut| is enough to indicate CP violation. Furthermore, 
the sinusoidal dependence implies as well a degeneracy be- 
tween Sfn — > 180° — 5/i t, as reflected in the figure. 

The right panel in Fig. Q] depicts the 3(7 sensitivities to 
|*7ut| (solid line) and S^r (dotted line), while the present 
bound from t — » fly is also shown (dashed line). The 
poorest sensitivity to around 10~ 3 , is found in the 

vicinity of S^ t = and <5/i T = 180°, where the CP-odd 
interference term vanishes and the bound is placed through 
the subleading | J7/it [ 2 term. The latter is also present at zero 
distance and its effects were already considered in Ref. []]], 
obtaining a bound of similar magnitude. The sensitivity to 
Itj^xI peaks around \ r\fi-z \ ~ 4 - 10~ 4 for <5/i T ~ ±90°, where 
sin5/i T is maximum. That is, for non-trivial values of <5/i T 
not only CP-violation could be discovered, but values of 
I i7/i t I an order of magnitude smaller could be probed. 

NON-UNITARITY VS 
NON-STANDARD INTERACTIONS 

Non-standard neutrino interactions (NSI) are usually intro- 
duced through the addition of effective four-fermion oper- 
ators to the SM Lagrangian 0,11, Q], made out of the SM 
fields and invariant under the SM gauge group. They can 
affect the production or detection processes or modify the 
matter effects in the propagation, depending on the oper- 
ator or combination of operators considered. For instance, 
in Ref. |4] new matter effects stemming from the operator 
{v S a' 'YuPlV^ (fYiif) were studied. In Ref. d, another 
kind of four-fermion operators were introduced which af- 
fect the production and detection processes at a Neutrino 
Factory, but do not correct matter effects. In that case, 
defining V™ = U a j V,-, with U being the PMNS matrix, the 




FIGURE 1. Left: 3d contours for two input values of \t]p.t\ and dfit represented by the stars. Right: the solid line 
represents the 3(7 sensitivity to |tJut| as a function of 8^, the dotted line the 3d sensitivity to S^t and the dotted-dashed 
line represents the present bound from T — » [iy. 



effective production and detection states are given by: 

\ v p> = (i+e*n e p\ v f>=(i+ e *n e pu*p i \ Vi > 

|v£> = (l + e* rf ) M/i |v™>=(l+ e w ) M/i ^,.|v / >, 

where are the coefficients of the new operators, up to 
normalization factors. These expressions are very similar 
to our parameterization of the effects of a non-unitary mix- 
ing matrix, Eq. ^3). In fact, the latter can also be encoded 
in terms of effective four-fermion operators, after integrat- 
ing out the W and Z bosons. The difference is that, in the 
case of a non-unitary mixing matrix, the coefficients of the 
different effective operators induced by it and contributing 
to production, detection and matter effects are not indepen- 
dent but related. It follows from Eq. that r) a p = T]p a . 
Which would mean = e^* , a constraint usually not re- 
quired when introducing NSI. When such equality holds, 
the oscillation physics induced by NSI is equivalent in 
vacuum to that stemming from non-unitarity in the MUV 
schemfl Furthermore, even if no such relations among 
and e^p are assumed, the order of magnitude of the bounds 
obtained in the previous chapters should apply as well to 
NSI, barring very fine tuned cancellations. The new CP sig- 
nals analyzed in this work are, therefore, also probes of the 
phases of NSI. 

CONCLUSIONS 

An asymmetry between the strength of — > V T oscilla- 
tions versus that for — » v T has been shown to be a beau- 

2 In Ref . 1 1 ] it was shown that matter effects are also modified in 
the presence of a non-unitarity mixing matrix, nevertheless, as we 
have studied a setup such that matter effects - standard and new 
ones - are negligible, the conclusion is the same even in matter. 



tiful and excellent probe of new physics, when measured at 
short-baselines (~ 100 km.) using a Neutrino Factory beam 
of energy ff{2QGeV). Non-trivial values of the new phases 
can lead not only to the discovery of CP-violation associ- 
ated to the new physics, but also allow to probe the absolute 
value of the moduli down to 10~ 4 . 

Our analyzes for future sensitivities to non-unitarity, as 
well as the new signals of CP- violation explored here, also 
apply to the physics of "non-standard neutrino interac- 
tions", except for extremely fine-tuned cancellations. 
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